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Abstract. Based on the conventional sequential-accretion paradigm, we have proposed that, 
during the migration of first-born gas giants outside the orbits of planetary embryos, super Earth 
planets will form inside the 2:1 resonance location by sweeping of mean motion resonances (Zhou 
et al. 2005). In this paper, we study the subsequent evolution of a super Earth (mi) under the 
effects of tidal dissipation and perturbation from a first-born gas giant (7712) in an outside orbit. 
Secular perturbation and mean motion resonances (especially 2 : 1 and 5 : 2 resonances) between 
mi and m2 excite the eccentricity of mi, which causes the migration of mi and results in a hot 
super Earth. The calculated final location of the hot super Earth is independent of the tidal 
energy dissipation factor Q' . The study of migration history of a Hot Super Earth is useful to 
reveal its Q' value and to predict its final location in the presence of one or more hot gas giants. 
When this investigation is applied to the GJ876 system, it correctly reproduces the observed 
location of GJ876d around 0.02AU. 
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1. Introduction 

The search for habitable planets is an essential step in the quests to unravel the origin of 
the Solar System and find life elsewhere. To date, more than 250 exoplanets are detected 
mainly by radial velocity survey of nearby solar-type star^. In the database, there are 
17 planets with mass less than 25 Earth mass (Mq), and among them 8 planets have 
orbits with period < 10 days. They are GJ876d, HD69830b, GJ674b, HD160691d, 55Cnc 
e, G1581b, HD219828b and GJ436b. We call them 'hot super Earths'. 

According to the conventional core-accretion scenario of planet formation, planets form 
in a protoplanetary disk around the host protostar. Through the sedimentation of dust, 
cohesive collisions of planetesimals, many embryos will form by accreting and clearing 
the planetesimals in their feed zone (a band centered on the embryo with a width of ~ 10 
Hill radius) and result in dynamically isolated bodies. In a disk with several (/d) times 
of minimum mass solar nebular, the isolation mass is (Zhou et al. 2007), 

A/iso = 0.51 X 10-^M(svkfi\ (1.1) 

where 

log(fciso) = Vb'^ + 0.61c - b, Q 2) 

6 = 2.8 -h 0.33 log 77, ^ ■ ^ 

c = 3.6 + 0.67 log ?7 + log Tdop, 



f http://vo.obspm.fr/exoplanetes/, http://exoplanets.org/[ 
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/ico — 1 for embryos inside the ice Hne and /ice = 4.2 outside that, Tdcp is the timescale 
of depletion of gas disk. 

According to equation (|l.ip . the isolation mass inside the ice line (o « 2.7AU in solar 
system) is too small to become a super Earth. Unless their nascent protostellar disks 
are highly compact, the observed super Earths are unlikely formed in situ. Some extra 
mechanisms are required to account for the excitation of eccentricity and the merge of 
isolated embryos into super Earths. In Zhou et al. (2005), we have proposed two mecha- 
nisms that may lead to the excitation of eccentricities of embryos: (1) During the type-II 
migration of a first-born gas giant planet outside the orbits of embryos, the locations of 
its mean motion resonances (mainly 2:1 resonance) sweep through the embryos region; 
(2) During the dispersal of the gas disk, the location of secular resonance between the 
gas giant and embryos sweeps through the inner orbits. Additional mechanisms have also 
been discussed by Raymond et al. (2007). 

In this paper, we suppose a super Earth has formed through one of the above mecha- 
nisms, and study the subsequent evolution after the gas disk was depleted and the gas 
giant has stopped its migration. First we briefly review the secular evolution of two plan- 
ets under tidal dissipation. Then we show some numerical results in section 3. Conclusions 
are presented in section 4. 



2. Secular dynamics under tidal dissipation 

2.1. Tidal perturbation timescale 

We adopt a two-planet system as a model. Suppose two planets with mass mi{i = 1,2) 
(in the order form inner to outer) moving around a star with mass in the same orbital 
plane. Let mi be an Earth-like planet, and TO2 a gas giant, Si,fli, the radius, 

spin rate (with spin axis perpendicular to the orbital plane), semi major axis, distance 
from the star of planet i (i = 1, 2respectively) , . The acceleration to the relative motion 
of mi caused by the tidal interaction between the star and planet mi has the form of 
(Mignard 1979, Mardling & Lin 2002) 




where f, are the unit vector of radial and transversal direction of the orbital plane, 
Vi = Virf + and ni are the Kepler velocity and mean motion of planet i(i=l,2), 
respectively, Q'^ and Q'^ are the effective tidal dissipation factor of the star and planet i 
defined as Q' = 3Q/{2kL), where — tan(2e) is the effective dissipation function, e is 
the tidal lag angle (Goldreich & Soter 1966), fc^ is the Love number or twice the apsidal 
Constance for gaseous planets(e.g., Mardling & Lin 2002), and 

^=&){—ri^)' (2.2) 

^i mi 

is the ratio of tidal dissipation in the planet to that in the star. If A 3> 1, tidal dissipation 
in the planet dominates the evolution. 

The values Q'^ inferred form the observation of circularization period in various stellar 
clusters are ^ 1.5 x 10^ for young stars with age less than O.lGyr, and ^ 10^ for mature 
stars (Terquem et al. 1998,Dobbs-Dixon et al. 2004). The Q' value for Jupiter inferred 
form lo's orbit evolution ranges from 5 X 10-* to 2 X 10^ (Yoder & Peale 1981). And for 
Earth, Q'^ « 60 (Yoder 1995). Thus for a gas giant planet with Jupiter mass, suppose 
Q'* ~ Q'j = 10^, from Ea. (|2.2p . A ^ 10, while for a terrestrial planet with Earth mass. 
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A ~ 10'*. So in the case of tidal interaction between an Earth-like planet and a star, tidal 
dissipation in the planet dominates the evolution, thus we neglect the contribution of 
tide in star in the following study. 

Take mi as an example. Under the perturbation of tidal effect, the averaged equations 
(over a period of orbital motion) governing the evolution of planet mi are. 



< ai >tidc= ^SoiTjij, 

< ei >tidc 



tide 

-1 
tid( 



/3(ei)-Ti(^)/4(ei) 

< tUi >tidc = < Al >tide= 0. 



/i(ei)-(^)/2(ei) 



(2.3) 



where vji, Ai are the longitude of perihelion and mean longitude of the orbit of nii (with 
volume density pi), respectively, and 



Ttide 



'^'^(^)(^)^=^2.4xl0Vr( 



Oi 



Pi 



63ni 5*1 
Functions used are 



13 . 3 mi 2 

) 2 ( ) 2 ( ) 3 ( - 

O.IAU^ Vq^ ^TOe ^3gcm- 



2\15/2 



hie) = (1 + fe' + ife^ + + f e«)/(l - e^) 

/,(e) = (l + fe2 + fe4 + Ae6)/(l-e2)6, 

/3(e) = (1 + f + f + l^e^)/{l - e2)W2, 
/4(e) = (l + |e2 + ie4)/(l- 6^)5, 
/5(e) = (l + 3e2 + |e4)/(i_e2)9/2, 
/6(e) = (1 + J^e^ + f|e4 + fe^ + ^e^ + ^e^')/{l 
Me) = (1 + + 8e4 + ffe^ + |f|e« + Tffcei°)/(1 



14'' ' '224^ ' 448 

The evolution of spin rate fii is subjected to, 

m*mi 



+ 3e^ + |e4), 
+ 3e2 + |e4). 



-ri X Fi, 



tide 



where h 



+ mi 

IniiSf is the inertial momentum of mi. The averaged change rate is 



yr- 

(2.4) 



(2.5) 



(2.6) 



< ^1 >tide" 2'^tid°(;5r)^ 



/2(ei)-(^)/5(ei) 
ni 



(2.7) 



A stable equilibrium configuration occurs at 



O _ /2(ei) 

"Leg — "TT 



(2.8) 



Since the timescale to reach the equilibrium state (~ Ttidc(*S'i/ai)^) is several orders less 
than the tidal circularization timescale, we suppose such a state is reached. Substitute 
Eq. (|2.8p into (|2.3p . we derive the timescales of tidal evolution of mi. 



-tide 



Ol 



(1 



22)15/2 



2e2/6(ei) 



'n;idc 1 '^c — tide 



ei 
ei 



(1 - e2)i3/2 
/7(ei) 



Ttide 



(2.9) 



Note that, Ta^ude > Te-ude whcu Bi « 0. Howcver, when ei w 1, Ta-ude and T^-ude 
could be very small, and Ta^ude < Te^ude as long as e > 0.63425.... 

Due to the huge difference of Q' between the Earth-like planet mi and the gas giant 
m2, for our later investigation of tidal evolution with ai < 0.63a2, we neglect the tidal 
effect in planet m2. 
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2.2. Secular evolution in the case of ei ^ 62 

When ei <^ 62, the secular evolution of mi , 7712 under tidal dissipation and general 
relativity effect can be approximated by the following equations (Mardling 2006): 



ei — — Woe2 sin?7 — WtCi, 
62 = WcCi sin?/, 
fj^Wq- Wo(f^)cos?7, 

where rj — vui ~ -cu2, a ~ ai/a2, l3 = — e|, and 



(2.10) 



(2.11) 



= fni(l2^)a3;3-3[l 



with 7 = 4(niai/c)^(m*/m2)a'^, the ratio of general relativity to quadruple contribution 
of 77. According to these equations, the secular evolution of ei and 62 mainly passes three 
stages: 

(1) After a short time oscillation, the evolution of ei and 77 reaches a state of librating 
around a quasi-equilibrium configuration with ei — e'j* and i] = 2mr or (2n-\- l)7r , where 
(Mardling 2006) 



62- 



w„ 



5/4ae2 



/32|l-^(mi/m2)/3-i+7/?3r 



(2.12) 



(2) As J] librates and ei evolves to ei = gradually, ai is damped according to Eq. 
thus mi migrates inward efficiently. 

(3) Finally 62 is damped on a timescale Tc ^ Te-ude- During this timescale, the orbit 
of mi is almost circularized, and migration of mi is effectively stopped at a location aif. 

The location of aif is what we want to find. However, due to the presence of resonant 
motion, the evolution of the two-planet system in real situation is more complicated, as 
we will show below. 




-1.4 -1.2 

log a,„ (AU) 

Figure 1. Maximum eccentricity (red solid line ) of mi in an initial circular orbit of semi-major 
axis flio excited by m2 (black circle, with 020 = O.IAU, 620 = 0.1). The black dotted lines 
with labels 2.0,4.0,... denote the timescale (log(ra_ti£je/years)) of mi from Eq. (|2.9|l at the 
specific location of (aio,ei) with Q'l — 100. The blue dashed line is obtained by two times of 
the equilibrium values defined by Eq. (|2.12|l . 
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3. Numerical simulations 

We study the migration of an Earth-like planet under the tidal and mutual planetary 
perturbations with general three-body model. The system consists of a solar mass host 
star (m* = IM©), an Earth- like planet with mass (toi = SAfg), and a Jupiter mass gas 
giant (m2 — Mj, where Mj is Jupiter mass). Let mi be initially in a nearly circular 
orbit (eio = 10^"^), and m2 with initial elements 020 = 0.lAU,e2o = 0.1. To shorten the 
integration time, we let Q' = 0.02, as the migration timcscalc is proportional to Q' . 

By integration of the full equations of the general three-body system without tidal 
dissipation, we plot the maximum eccentricity (eimax) of mi excited by m2 in Figure 1. 
The corresponding tidal-damping timescale obtained from equation (|2.9p with (aio,ei) 
is also shown in the background of Fig.l. 

As we can see from Figure 1, the orbits of mi at most locations with aio < 0.16AU have 
7a-tido < lOGyr. However, most of those orbits in the Hill unstable region around m2 with 
half width A = (e2o + 2V3/i)a2 « 0.034AU (where h = [7712/(3771*)]^/^) wiU be scattered 
to far away in our coplanar model. According to Zhou et al. 2005, embryos formed inside 
the location of the 2 : 1 resonance(a2:i « 0.063AU) with 7^2 are dynamically stable, so 
we focus on the evolution of orbits with initial semi major axis oio ^ 0.063AU. 

If planet mi is initially located in lower order mean motion resonances with 7712, its 
eccentricity will be excited, thus a fast inward migration of mi is induced, according to 
Eq. (|2.9p . Fig. 2 shows the evolution of two orbits either from 2 : 1 resonance location 
(oio = 0.063AU) or from non-resonance location (aio = 0.040AU). During the subsequent 
passage through 5 : 2 resonance, the amplitude of eccentricity excitation is relatively 
large. Recall that, according to Eq. (j2.9p . the a— damping timescale is much smaller than 
that of e-damping at high eccentricity. Thus a fast migration occurs until ei decrease to 
a small value ~ 0.01 (Novak et al. 2003). Then a slow migration linked with the secular 
dynamics occurs, with 77 = zui — wi librating along an equilibrium value(see Fig. 2b). 
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Figure 2. Evolution of orbits with aio = 0.063 (orbit A) and aio = 0.040 (orbit B). (a) Evolution 
track in ai — ei plane. The dashed line shows the width of 2 : 1 resonance obtained from the 
circular restricted-three-body problem. The 5 : 2 indicates the 5 : 2 resonance location at the 
place that the eccentricity of orbits A jumps up. The arrows indicate the evolution directions, 
(b) Evolution of — vj^) (in radian). 
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The migration induced by resonant eccentricity-excitation is different from that excited 
by mutual secular perturbation. When we check the evolution of 0^,6^,(1 — 1, 2) during 
the passage of 5 : 2 resonance, we find that 01,61,62 have dramatic decrease after the 
crossing the 5 : 2 resonance (Fig. 3). The decrease of 62 causes the different final states 
(i.e., the final state of orbit A and B) in (ai,ei) plane of Fig. 2. The final location is 
around 0.018 ~ 0.025AU, depending on the different evolutionary routines. 

In order to show that the above track correctly reproduces the observed location of 
extrasolar planets, we applied this study to the GJ876 system. GJ876 is a M dwarf 
star located 4.72 pc away from us in the solar neighborhood. To date, two gas giant 
planets, GJ876b and GJ876c, were observed to be located on orbits with period around 
30 days and 60 days, an example of 2:1 mean motion resonance, and a hot planet GJ876d 
with mass around 5.7Mq in an orbit with period 1.94 days(a = 0.0208AU, e=0). We 
numerically simulate the evolution of an Earth-like planet inside a gas giant located in 
the present orbit (020 = 0.13AU, 620 0.2243). Figure 4 shows the evolution track of 
mi. The final location is around 0.02 AU according to the simulation, which is almost 
independent of the initial location of mi. 



4. Summary and discussions 

Many super Earths are observed to be located inside the orbits of gas giants. These 
super Earth and gas giant pairs may be a natural consequence of planet formation and 
migration. Embryos formed prior and interior to the gas giants are induced to migrate, 
collide, and evolve into close-in Super Earths(Zhou et al. 2005). In this report, we have 
shown that, the migration of super Earths under tidal dissipation and the perturbation 
from gas giants is mainly along the secular evolution paths. Although resonances be- 




Figure 3. Evolution of orbital elements for the two orbits shown in Fig. 2 with initial elements 
aio = 0.040, 0.063, respectively. The evolution of orbit with aio = 0.040 fits with the secular 
evolution described in section 2, but for the orbit with aio ~ 0.063, the presence of the resonance 
leads to a dramatic increase (decrease) of 61(62) at time t ~ 8x lO'* year, when the orbit crosses 
the 5:2 resonance in Fig. 2. 
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Figure 4. Evolution track of orbits witli aio = 0.07074 (tlie location of 5 : 2 resonance, orbit 
A) and aio = 0.060 (orbit B). The arrows indicate the evolution directions. 

tween super Earths and gas giants may excite the eccentricity and speed the migration 
timescale, the final evolution path can be well determined. 

According to the investigation of this paper, we find that the study of the evolution 
path provides useful information in the following ways: (i) the migration path shows 
the evolution history, especially the evolution of ei, 62, ai (Figs. 2, 4). (ii) Comparing the 
observed location of the planet in the path, we can deduce the range of Q' values for the 
hot Super Earths. We will investigate in more details on these topic in the future. 
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